program bomb
implicit none
real	::	h=1000.0,v=0.0,t=0.0,dt=0.02

open(1,file='height.dat')
do while(h>0)
	v=v+dt*dv(v)
	h=h-dt*(v)
	t=t+dt
	write(1,*)	t,h
end do

contains
	function	dv(v)
	real,	intent(in)	::	v
	real				::	dv
	dv=9.8-0.225*v*v/90
	end function
	
	function	dx(v)
	real,	intent(in)	::	v
	real				::	dx
	dx=v
	end function


end program
